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ABSTRACT 

Aims. We present mid-IR spectro-interferometry of the Seyfert type 1 nucleus of NGC 3783. The dusty circumnuclear environment 
is spatially resolved and the wavelength dependence of the compact emission is discussed. 

Methods. The observations were carried out with the MIDI instrument at the Very Large Telescope Interferometer in the ,/V-band. 
Spectra and visibilities were derived with a spectral resolution of A/AA ~ 30 in the wavelength range from 8 to 13^m. For the 
interpretation we developed a simple dusty disk model with small and variable covering factor. 

Results. At baselines of 65 and 69 m, visibilities in the range of 0.4 to 0.7 were measured. The ,/V-band spectra show a monotonic 
increase of the measured flux with wavelength with no apparent silicate feature around 10 pm. We find that the mid-IR emission from 
the nucleus can be reproduced by an extended dust disk or torus with a small covering factor of the radiating dust clouds. 
Conclusions. Our mid-IR observations of NGC 3783 are consistent with a clumpy circumnuclear dust environment. The interpretation 
in terms of a dusty torus with low covering factor supports a clumpy version of the unified scheme for AGN. The inferred sizes and 
luminosities are in good agreement with dust reverberation sizes and bolometric luminosities from optical and X-ray observations. 
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1. Introduction 

Unified scenarios of AGN suggest that the central continuum 
source is surrounded by a dusty molecular torus. This torus is 
thought to be responsible for a spect-angle-depe ndent obscura- 
tion of the central source (e.g.. lAntonuccil l 1993). Mid-infrared 
(MIR) inte rfero metry of the Seyfert 2 nucleus in NGC 1068 
dJaffe et alj|2004t iPoncelet et alj|2006l) and the Circinus galaxy 
dTristram et al.l2007h provided supporting evidence for the dusty 
torus model. The MIR emission of these nuclei was resolved by 
MIR interferometry with the MIDI instrument in both cases. The 
resolved structures are belived to be the IR-bright innermost re- 
gion of the torus. In Seyfert type 1 sources the situation is not as 
clear as NIR interferome try data obtained with the Keck interfer- 
ometer on NGC 4151 dSwain et alj|2003l) have been interpreted 
as accretion disk emission. 

At a redshift of z = 0.0097 dTheureau et al.1 11998L D = 
44Mpc), the galaxy NGC 3783 harbors a bright Seyfert 1 
core with a black hole of M B h = 3.0 ± 0-5 x 1O 7 M derived 
from reverberation mapping dOnken et alJl2004h . NGC 3783 has 
been the target of ex tensive ultraviolet spectral monitoring (e.g., 
lReichertetal]ll994l). A narrow iron K line is seen in XMM- 
Newton spectra dReeves et alj 120041) consistent with an origin 
in Compton-thick matter such as a molecular and dusty torus. 
The optical spectrum shows blueward asymmetries in forbid- 
den lines and broa d lines wi th up to 3800 km/s in the broad- 
est components dEvansI U988I) . T he extended narrow line re- 
gion of NGC 3783 seen with HST dSchmitt et alj|2003l) appears 
halo-like in [Om] w ith a size of 1"9 (400 pc). The NIR spec- 
trum of the nucleus dReunanen et alj|2003l) features a broad Bry 
line and three coronal lines are detected. M I R spe ctra within a 
beam of 4'.'2 were presented bv lRoche et al.l d!991l) and showed 
an almost featureless spectrum with a V-band luminosity of 
~ 3 x 10 9 L G in this aperture. VISIR photometry presented by 



iHorst et alj d2006l) show that NGC 3783 falls onto the MIR vs. 
hard X-ray correlation for AGN with a rising flux of 0.66 to 
0.72 Jy across the V-band. The UV luminosity is estimated from 
Luv ~ 6 AL A (5 100 A) to b e 1 - 2 x 10 37 W. T his is consitent 
with the Lboi - Lx-relation ( Marconi et alj |2004) which predicts 
a bolometric luminosity of 3 x 10 37 W for NGC 3783. The accre- 
tion disk around the central black hole is therefore radiating at 
~ 10% of the Eddington limit. All these characteristics of a clas- 
sical type 1 core of suggest the presence of a dusty nuclear torus 
seen almost face-on. Direct evidence of this torus can only be ob- 
tained through high spatial resolution interferometry of the hot 
dust emission in the near- and mid-IR. In this paper we present 
the first spectro-interferometric observations of NGC 3783 with 
MIDI/VLTI. In Sec. 2 we describe the observations and data re- 
duction. In Sec. 3 we develop a simple model for the circumnu- 
clear dust distribution and discuss the implications for dusty tori 
in the unification scheme of AGN. Sec. 4 gives a short summary. 



2. Observations 

NGC 3783 was observed with MIDI in the two nights 27 and 
30 of May 2005 (Program 75.B-0697). A prism with a spectral 
resolution of A/AA = 30 at 10 fim was used to obtain spectrally 
dispersed fringes between 8 and 1 3 [i m. A detailed descriptio n 
of the observing procedure is given in Przvgoddaet al. (2003). 
We use the data reduction package EWS developed at Leiden 
Observatory (combined MIA+EWS, ver.l.5.2fl The EWS soft- 
ware first corrects for instrumental and atm ospheric p ath differ- 
ences before adding the fringes coherently (Jaffe 2004). We de- 
rive the interferometer transfer function at each wavelength (i.e., 



1 Available at http://www.mpia-hd.mpg.de/MIDISOFT/ and 
http://www.strw.leidenuniv.nl/nevec/MIDI/index.html 
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Table 1. List of calibrators used in this work, together with 12 
yum fluxes (-F12), and uniform-disk diameters ( duo) of the MIDI 
observations. 



Calibrator 


Fu (Jy) 


d VD (mas) 


HD 100407 


11.6 


2.41 


HD 116870 


10.0 


2.57 


HD 123139 


54.4 


5.25 


HD 169916 


22.4 


3.90 


HD 160668 


28.8 


2.30 


HD 152885 


12.6 


2.88 


HD 178345 


10.6 


2.50 
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Fig. 1. Visibilities vs. wavelength of NGC 3783 in the 8 - 13 /mi 
range observed with MIDI on 30 May 2005. Visibilities from 
dataset 1 (diamonds) show smaller visibilities than dataset 2 
(crosses) for projected baselines of B — 68.6m at P.A.= 115° 
and B = 64.9 m at P.A.= 120°, respectively. The gray area be- 
tween A = 9.4 and 9.8 /mi is affected by the atmospheric ozone 
band, which could not be properly calibrated. The model vis- 
ibilities A-l (blue, dash-dotted) A-2 (green, dash-dotted), B-l 
(blue, solid) and B-2( black,solid) are shown for comparision. 
The models A-2 (green) and B-2 (black), modified for the longer 
baseline of B — 68.6 m, are also plotted as dotted curves. This 
demostrates that the visibilty offset between No. 1 and 2 is not 
explained by the changing baseline when the visibility fo each 
dataset utilizes the associated photometry shown in Fig. [2] 



at each spectral channel) between 8 and 13 /mi by observing cal- 
ibrators whose angular diameters are known (see TableQ]). 

In order to calibrate the visibility of NGC 3783 at a given 
wavelength, we use the mean of the transfer function values 
derived from all the calibrators observed on the same night as 
NGC 3783. The errors of the calibrated visibilities are esti- 
mated from errors on the raw visibilities and the standard de- 
viation of the transfer function values at each wavelength (added 
in quadrature). The relative errors of the calibrated visibilities 
are typically ~ 10%. Fig. [TJ shows the calibrated visibilities of 
NGC 3783 from two datasets (No. 1&2 in the following) ob- 
tained close in time in the night 30 May 2005. The observed 
TV-band visibilities are in the range of 0.4 to 0.7. The wavelength 
dependence is weak within the errors from 8 to 1 3 /mi for both 
datasets. The projected baselines are B = 68.6 m and 64.9 m 
for the dataset 1 and 2, respectively. We include in Fig. [TJ also 
the visibilities of two different models described in Sec. 13. 21 be- 
low. The wavelength region around 9.6 /mi (marked as a gray 
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Fig. 2. MIDI TV-band spectra of the nucleus of NGC 3783 in Jy 
vs. wavelength. The spectra have been calibrated and averaged 
against the reference stars of the night 30 May 2005. Dataset 
1 and 2 are shown with diamonds (1) and triangles (2), respec- 
tively. The model spectra for B-l (blue) and B-2 (black) of a 
clumpy and dusty torus corresponding to Fig. 1 and described 
in Sec. 13. 21 are overplotted. The weak or absent silicate feature 
around A = 9.7 /mi is a consequence of the large optical depth in 
the models. 

area in Figs.QJandO is difficult to calibrate for faint targets with 
F v ss 1 Jy due to a strong atmospheric ozone band. This may ex- 
plain the apparently larger visibilities in that wavelength region. 
In addition, we extracted absolute TV-band spectra of NGC 3783 
from MIDI single channel photometry Interferometric calibrators 
observed in the same night as NGC 3783 and at similar air- 
masses serve as spectrophotometric standard stars to calibrate 
the spectra (see Table [TJ). The cali brated spectra of these stars 
are taken from ICohen et alj d!999l) . The two calibrated spectra 
of NGC 3783 are shown in Fig. [2] We estimate the errors of 
the spectra to be 10-20%. The flux spectrum rises slowly with 
wavelength over the whole 8-13 /mi ran ge. It is important t o 
note that neither the flux, in agreement with lRoche et alj (1199 11) . 
nor the visibility spectrum shows a prominent silicate feature at 
~ 10 /mi. 

Due to variable thermal emission from the sky and from in- 
side the interferometer, the uncorrected photometry measure- 
ments are subject to larger errors than the correlated flux mea- 
surement. In an alternative approach we combined the spectra 
from both datasets to create a common total flux spectrum shown 
in Fig. [3] We used this to derive corrected visibilities for datasets 
1 and 2 also included in Fig. [3] A comparision with model B- 
3 shows that the remaining difference between the two datasets 
can be attributed to the difference in the baseline. 

3. Discussion 

3. 1 . General properties 

The first and most notable feature of the visibilities in Fig.Q] 
and the spectra in Fig. [2] is the absence of the 9.7 /mi absorp- 
tion or emission feature of silicate dust expected to be present in 
the torus. This absence may be due to missing silicates in AGN 
dust, because of the lower sublimation temperature Tsi ~ 1000K 
when compared with graphite grains Tq ~ 1500 K. But the mid- 
IR emission is expected to come from T ~ 300 . . . 800 K dust, 
where both kinds of dust grains should be present. Alternatively, 
the dust grain size distribution maybe dominated by large grains, 
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Fig. 3. Combined MIDI /V-band spectrum (top) and visibilities 
(bottom) of NGC 3783. We used an arithmetic mean with errors 
for the two spectra from Fig. [2] and compare it with model B- 
3 (blue solid curve). The visibilities for datasets 1 (diamonds) 
and 2 (crosses) have been derived with the combined spectrum 
shown at the top. The visibilities for model B-3 for baseline 
B = 68.6 m (solid) and 64.9 m (dashed) are included. The mea- 
surements deviate substantially from the model in the 8.3 to 
9.3 /mi region. 



for which the silicate feature is less prominent. We have no di- 
rect evidence for this in NGC 3783, except for dust-reverberation 
arguments (see below). The third possibility, and our choice for 
the models below in Sec. 13. 21 is a radiative transfer effect: For 
moderate optical depth in the 8-13 fim range with t ~ 0.5 and 
small temperature variations along the line of sight through the 
dusty torus, the emission spectrum approaches a black body and 
the feature is suppressed. 

Secondly, the flux spectrum F v in Figs.|2]and[3]shows a linear 

rise with wavelength, in contrast to most radiative transfer simu- 

j ii 1 

lation of face-on dusty AGN tori (Honig et al. 2006, e.g.,). This 

may indicate the presence of additional, spatially distributed 
heating sources, which can increase the amount of dust at in- 
termediate temperatures of 150 to 300 K. Another possibility is 
a smaller covering factor of hot dust compared to cooler dust in 
a disk heated by the central engine. This idea is employed in the 
models below. 

In addition, we find almost constant visibilities between 8- 
13/im suggesting a linear rise of the source size with wave- 
length. The derived visibilities are smaller than expected for a 
compact AGN dust torus, for which the NIR and MIR emission 
is dominated by a small region close to the dust sublimation ra- 



dius. The diameter of the inner rim of the torus is expected (see 
Eq.|2]below) to be 2.1 mas corresponding to V close to unity. A 
small dust covering factor in the models below can explain the 
data and imply a large extension of the MIR emission region. 

Information on the torus inclination can be infered from 
spectropolarimetry in the optical 400 to 500 nm region, which 
show s a relatively larg e polarisation fraction of ~ 1 % at R A. ~ 
130° dLira et al.ll2.007l and Kishimoto priv.comm.). The axis of 
the optical scattering region is therefore inclined with respect to 
our line of sight suggesting the minor axis of an elongated MIR 
disk to be close to the P.A.~ 120° of our observations. The RA. 
of the interferometric baselines for the datasets 1 & 2 are very 
similar, so that we cannot constrain any non-circularity of the 
MIR emission from our data. 



3.2. Modeling visibilities and spectra with a dusty disk 

For the interpretation of the data presented a bove we employ a 
simpl i fied version of the clumpy torus model dBeckert & Duschll 
2004; iHonig et af] l2006). A torus seen almost face-on (z « 0°) 
appears as a disk composed of dusty clouds. The inclination has 
almost no observable effect on spectra and colors as long as i 
is smaller than the half opening angle of the torus. We assume 
single temperature clouds, where the cloud temperature depends 
only on the radius in the disk. F or simplicity and c omparison 
with other models we use a MRN (Mat his et al.ll 19771) grain size 
distribution between 0.005 fim and 0.25 fim. The inner radius of 
the disk or torus is at R[ n , where we assume a dust temperature 

of Ti n 1500 K, which is s uggested by NIR colors of type 1 

AGN dKishimoto et al.ll2007h . The local intensity / v (r) emitted 
by clouds is a modified black body B(v, T) at temperature T = 
Tin x (f/Rin) ' 1 with optical depth r y , so that 



I v (r) = B(v,T(r))(l-e T ') 



(1) 



We are interested in lines of sight through the disk parallel to 
the disk symmetry axis. For simplicity, we use a constant t v 
independent of the radius. The frequency dependence of r v is 
derived from the MRN distribution and a mixture of 50% sili- 
cates and 50% graphite grains wit h the optical properties taken 
from Weingartner & Draine (2001, and references thereinfl For 
comparison we fix t v in the optical Ty at A = 560 nm. For 
normal black body emission in radiative equilibrium with the 
central illuminating UV/optical radiation we expect an exponent 
77 = 0.5. This leads to an almost constant flux spectrum F v when 
integrated over the disk surface. We employ two modifications 
(models A and B) to fit the total flux spectrum in Fig. [2] (dash- 
dotted lines). Model A-l & A-2 for dataset 1 & 2 assumes a dif- 
ferent exponent 77 = 0.35 (see Table|2]i. This is motiva ted by the 
equilibrium temperature of graphite grains (see e.g., iBarvainisI 
119871; iKishimoto et al.ll2007l) 



1500 K 



10 37 W 



0.179 



R 



0.13pc 



-0.357 



0.05 fim 



-0.179 



(2) 



where the grain size a = (a 3 ^ / (a 2 ^ is an average size for the 
grain size distribution. Such a model can explain the spectra in 
Fig|2]but not the visibilities as shown in Fig.Q] We can accom- 
modate the suggested large size of the dust emitting region by 
introducing a covering factor f(r), for which we assume a ra- 
dial dependence f(r) = fo(r/R m )^. For models A, we use £, — 
and /o can be adjusted to get the size of the source. The idea of 



2 Files with tabulated optical constants have been taken from 
http://www.astro.princeton.edu/~draine/dust/dust.diel.html 
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Table 2. Dusty disk/torus model results. The covering factor /o 
and the exponents £ and rj for the covering factor and dust tem- 
perature are described in Sec. 13. 21 t> is the optical depth along a 
vertical path through clouds in the disk taken at 560 nm. 



Model 


TV 


/o 


( 


n 
















[pc] 


[10 37 W] 


A-l 


6 


0.04 





0.35 


0.064 


0.16 


A-2 


12 


0.06 





0.35 


0.04 


0.065 


B-l 


4 


4.0 10~ 3 


0.8 


0.5 


0.30 


3.7 


B-2 


8 


3.8 10" 3 


0.9 


0.5 


0.21 


1.7 


B-3 


8 


3.2 10~ 3 


0.9 


0.5 


0.23 


2.2 



Long term optical/IR photometry of NGC 3783 dGlassll2004 
and references therein) suggests a time delay of 190d between 
U- and L-bancQ. This implies a physical path difference of 
0.16pc or an angular scale of 0.75 mas, if the response is in 
the plane of the torus/disk. Given the uncertainties entering in 
Eq.|2]we get an agreement with the B-models within a factor 
of 2. For comparisio n, the NIR interferometry of NGC 4151 
by IS wain et al. I (120031) suggests a more compact emission region 
like an accretion disk or a ring at Ri n . A discussion of dust sub- 
limation radii, dust reverb eration time delay and luminosity de- 
pendence can be found in Kishim oto et al] (120071) . The inferred 
AGN luminosity from the B-models (Table© of L boi * 3 10 37 W 
or from the reverberation size via Eq.[2]are in nice agreement 
with the L b oi - Lx-relation dMarconi et a l. 2004). 



a low covering factor is a direct c onsequence of clumpy torus 
models (see e.g., Honig et al. 2006). In this scenario only clouds 
directly illuminated by the central source contribute significantly 
to the torus emission and due to the clumpy structure there is a 
non-zero probability for such clouds at all radii. 

The implied AGN luminosity for models A-l & A-2 are 
much smaller (see Tabel|2]i than derived from optical and X-ray 
observations. We therefore concentrate in the following on mod- 
els B-l & B-2 which provide a better fit to the spectra in Figj2] 
and the visibilities in Fig. Q] with a running covering factor. All 
models exhibit decreasing visibilities with wavelength in Fig.Q] 
and both models A and B give acceptable fits to the datasets 1 
& 2. The parameters of the models are summarized in Table|2] 
The inner torus radii of the models are determined by the total 
flux level and the covering factor in the combination /?;„/ JJq. 
The AGN luminosity (assuming isotropic emission) in Table|2]is 
derived from Eq.[2] The larger luminosities for B-l, B-2 and B- 
3 are a direct consequence of the smaller covering factor at the 
inner boundary R m . Note that the angular scale for R\ n = 0.1 pc 
is only 0.47 mas. The actual T = 300 K diameter is as large as 
70 mas for model B-l or 60 mas for A- 1 . 

The visibilities derived with the combined total flux spec- 
trum in Fig. [3] are much more consistent with each other. The 
difference between dataset 1 & 2 can be completely attributed 
to the difference in baseline length. This allows to fit the long 
wavelength side the NGC 3783 datasets with a single model 
B-3 which is very similiar to B-2 (see Table [2J. In the 8.3 to 
9.3 fim region the derived V{A) is substantially lower than sug- 
gested by all simple power-law models. NGC 3783 appears to 
be even larger than the models in this wavelength region. 

3.3. Implications for unified models 

Models o f dusty nuclear t ori suggested by the unified scheme 
of AGN (Anto nuccilll993t) r equire the dust to be conta ined in 
individual clouds or clumps dKrolik & Begelman] |1988l) . In re- 
cent years, models for the cloud size and cloud distribution 
dBeckert & Du schl 2004 ) and for the emission from a cl umpy 
dust distribution dNenkova et al .1120021: iHonig et al]|2006l) have 
been developed. The dyn amical model based on cloud collisions 
(Beckert & Duschll2 004) suggests that a geometrically thick and 
obscuring torus is expected for high mass accretion rates close 
to the Eddington limit of the central black hole. From the X-ray 
and optical/UV luminosity of NGC 3783 an Eddington ratio of 
10% has been inferred. Our simple model for the mid-IR emis- 
sion suggests that the 'torus' has a much smaller covering factor 
than expected for the collisional scenario. 



4. Summary 

We present MIDI/VLTI spectro-interferometry of the Seyfert 1 
nucleus of NGC 3783 and spatially resolve the object between 
8 and 13 fim. The observed visibilities suggest a significantly 
larger MIR emission region than the estimated dust sublimation 
size of 1 .4 mas. Spectra and visibilities are noticeable for the ap- 
parent absence of a silicate feature around 10/im. In addition the 
flux spectrum shows a linear increase with wavelength through- 
out the TV-band. Consistent results are obtained with the use of 
an average total flux spectrum (Fig. [3} which exemplifies that the 
incoherent flux measurement contributs the largest uncertainties. 
Our simple modeling reproduces simultaneously the observed 
visibilities and spectra by a dusty disk with small covering fac- 
tor seen face-on with a moderate optical depth. This model is in 
line with recent clumpy torus models, albeit the covering factor 
of clouds in NGC 3783, which contribute to the MIR emission, 
is very small and rises outwards in the disk or torus. 

Future near-IR and additional mid-IR interferometric mea- 
surements with different baseline length and covering a wider 
range of position angle are needed to confirm the present results 
and probe the validity of the simple models presented here. 

Acknowledgements. T.B. likes to thank Christian Hummel and Sebastien Morel 
at ESO/Paranal for help and advice during the observations. We like to thank 
Makoto Kishimoto for stimulating discussions and helpful comments. 



References 

Antonucci, R. 1993, ARA&A, 31, 473 

Barvainis, R. 1987, ApJ, 320, 537 

Beckert, T. & Duschl, W. J. 2004, A&A, 426, 445 

Cohen, M., Walker, R. G., Carter, B., et al. 1999, AJ, 1 17, 1864 

Evans, I. N. 1988, ApJS, 67, 373 

Glass, I. S. 1992, MNRAS, 256, 23P 

Glass, I. S. 2004, MNRAS, 350, 1049 

Honig, S. E, Beckert, T., Ohnaka, K., & Weigelt, G. 2006, A&A, 452, 459 
Horst, H, Smette, A., Gandhi, P., & Duschl, W. J. 2006, A&A, 457, L17 
Jaffe, W. 2004, SPIE Proc, 5491, 715 

Jaffe, W., Meisenheimer, K., Rottgering, H. J. A., et al. 2004, Nature, 429, 47 
Kishimoto, M., Honig, S. E, Beckert, T., & Weigelt, G. 2007, A&A, 476, 713 
Krolik, J. H. & Begelman, M. C. 1988, ApJ, 329, 702 

Lira, P., Kishimoto, M., Robinson, A., et al. 2007, in Astronomical Society of the 
Pacific Conference Series, Vol. 373, The Central Engine of Active Galactic 
Nuclei, ed. L. C. Ho & J.-W. Wang, 407-+ 
Marconi, A., Risaliti, G., Gilli, R., et al. 2004, MNRAS, 351, 169 
Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425 
Nenkova, M., Ivezic, Z., & Elitzur, M. 2002, ApJ, 570, L9 
Onken, C. A., Ferrarese, L., Merritt, D., et al. 2004, ApJ, 615, 645 
Poncelet, A., Perrin, G., & Sol, H. 2006, A&A, 450, 483 



In an earlier paper dGlassI 19921) a time delay of 85 d was suggested, 
which would imply a much smaller physical scale for the response. 



Thomas Beckert et al.: Probing the dusty environment of the Seyfert 1 nucleus in NGC 3783 with MIDI/VLTI interferometry 



5 



Przygodda, R, Chesneau, O., Graser, U., Leinert, C, & Morel, S. 2003, Ap&SS, 
286, 85 

Reeves, J. N., Nandra, K., George, I. M., et al. 2004, ApJ, 602, 648 
Reichert, G. A., Rodriguez-Pascual, P. M., Alloin, D., et al. 1994, ApJ, 425, 582 
Reunanen, J., Kotilainen, J. K., & Prieto, M. A. 2003, MNRAS, 343, 192 
Roche, P. F., Aitken, D. K., Smith, C. H., & Ward, M. J. 1991, MNRAS, 248, 
606 

Schmitt, H. R., Donley, J. L., Antonucci, R. R. J., Hutchings, J. B., & Kinney, 

A. L. 2003, ApJS, 148, 327 
Swain, M., Vasisht, G., Akeson, R., et al. 2003, ApJ, 596, L163 
Theureau, G., Bottinelli, L., Coudreau-Durand, N., et al. 1998, A&AS, 130, 333 
Tristram, K. R. W., Meisenheimer, K., Jaffe, W., et al. 2007, A&A, 474, 837 
Weingartner, J. C. & Draine, B. T. 2001, ApJ, 548, 296 



